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I . INTRODUCTION 


There  exist  phenomena,  such  as  solar  proton  events  and  nuclear 
explosions,  that  can  severely  disturb  the  total  charge  content  of  the 
atmosphere.  In  the  altitude  interval  30-80  km  these  disturbances  can 
affect  systems  over  a wide  range  of  the  electromagnetic  spectrum,  e.g., 
disruption  of  radar  and  communications  operations.  The  rate  at  which 
the  disturbed  atmosphere  returns  to  normal  depends  on  the  processes  for 
removing  the  excess  charge.  At  the  higher  altitudes  (70-80  km)  electron- 
ion  recombination  occurs  faster  than  negative  ion  formation  and  so 
charge  neutralization  is  dominated  by  the  electron-ion  dissociation 
recombination  process.  At  lower  altitudes  (<  ~ 60  km)  negative 
ion  formation  permits  ion-ion  mutual  neutralization  to  compete  with  and 
(still  lower)  to  dominate  electron-ion  dissociative  recombination. 

(The  precise  times  and  altitudes  at  which  this  dominance  occurs  is  a 
function  of  the  intensity  of  the  disturbance  as  well  as  the  rates  of 
the  recombination  processes.  In  this  report  we  shall  focus  on  the  ion- 
ion  mutual  neutralization  process  for  fixed  excitation  conditions  and  so 
reported  times  and  altitudes  would  have  to  be  adjusted  for  other  disturb- 
ing conditions.) 

For  ELF/VLF  communications,  Field  and  Dore1  have  found  that  the 
largest  uncertainties  in  these  attenuation  rates  (several  dB  per  megameter) 
can  be  attributed  to  uncertainties  in  the  positive  ion  - negative  ion 
recombination  rate  coefficient,  a..*  They  have  also  found  that  at  VLF 
uncertainties  in  the  ion-neutral  collision  frequency,  v. , can  cause 
uncertainties  in  attenuation  about  equal  to  that  found  £or  uncertainties 
in  a^. 


These  two  parameters  are  indirectly  related,  v.  is  a function  of 
the  ion  distribution.**  This  distribution  is  determined  in  part  by  the 
recombination  coefficient.  The  question  that  arises  is  whether  or  not 
the  operational  uncertainties  due  to  v.  and  to  can  be  treated  as 
separate,  independent  problems.  To  i estate  this  problem-  What  is  the 
change  in  the  ion  species  distribution  for  a given  change  in  a^? 


The  subscript  "i"  refers  to  "ion"  as  opposed  to  electron. 

* 2 -k 

For  one  ion  species  and  one  neutral  species  a.  <*  A p , where  A is 

the  sum  of  the  molecular  and  ionic  radii  and  p is  the  reduced  mass  of 

the  colliding  ion  and  molecule.  (Field,  private  communication  quoting 

Propagation  of  Electromagnetic  Waves  in  a Plasma  by  Ginzberg,  trans. 

by  Roger  and  Roger,  Gorden  and  Breach  Publishers,  N.Y.,  p.  97.) 

*E.  C.  Field  and  M.  N.  Dore,  private  communication,  1976. 
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Though  some  information  regarding  a.  has  been  determined  from 
analysis  of  natural  and  nuclear  disturbances,  ^uch  reliance  has  been 
placed  on  model  predictions.  The  AIRCHEM  code  has  been  used  at  the  BRL 
to  predict  the  response  of  the  ionosphere  to  disturbed  conditions  and  as 
such  it  provides  us  with  a convenient  tool  for  examining  the  sensitivity 
of  ionic  populations  to  variation  in  values  of  individual  two  body  ion-ion 
recombination  rate  coefficients,  a..  We  distinguish  between  an  effective 
ion-ion  recombination  coefficient,  ou  and  the  coefficient  foj  the  jth  + 
recombination,  a.,  a.  is  a weighted  sum  of  the  a.,  a = En.  n.  a./En.  n.  , 
where  n.  and  n.?  areHhe  positive  and  negative  ion  densities  of  the  J 
jth  recombination  reaction. 

This  report  describes:  our  current  state  of  knowledge  about  the 

values  of  a.  (section  II),  the  procedure  used  in  those  sensitivity  tests 
and  the  commuted  changes  in  the  ion  populations  for  changes  in  selected 
values  of  ou , under  fixed  disturbing  conditions  (section  III). 

II.  VALUES  FOR  THE  TWO  BODY  ION- ION  RECOMBINATION  COEFFICIENT,  a. 

3 4 

This  topic  has  been  reviewed  in  1972  and  1974,  and  will  be  dealt 
with  briefly  here.  Room  temperature  values  of  the  recombination  coeffi- 
cient for  fight  pairs  of  ions  of  atmospheric  interest  are  listed  in 
Table  I.  ” The  salient  feature  to  be  gleaned  fr<jim  Table  I is  tlje  wide 
range  of  disagreement  for  all  ion  pairs  except  O2  + NO^  and  02  + 0 . 

Values  differ  by  about  a factor  of  ten  in  two  cases. 

^E.  L.  Lortie,  M.  D.  Kregel  and  F.  E.  Niles,  "AIRCHEM:  A Computational 
Technique  for  Modeling  the  Chemistry  of  the  Atmosphere,"  BRL  Report 
No.  1913,  1976.  (AD  #A030157) 

3 

B.  H.  Mahan,  "Recombination  of  Gaseous  Ions,"  in  Advances  in  Chemical 
Physics  23,  1-40,  1973,  Ed,  by  I.  Prigcgine  and  S.  A.  Rice. 

4 

J.  T.  Moseley,  R.  E.  Olson  and  J.  R.  Peterson,  "Ion-Ion  Neutralization," 
Case  Studies  in  Atomic  Physics  5,  1-45,  1975. 

^M.  N.  Hirsh  and  P.  N.  Eisner,  "Laboratory  Measurements  of  Ion  Chemistry 
in  a Simulated  Disturbed  Ionosphere,"  Radio  Sci.  7_,  125-131,  1972. 

^D.  Smith  and  M.  J.  Church,  "Binary  Ion-Ion  Recombination  Coefficients 
Determined  in  a Flowing  Afterglow  Plasma,"  International  Journal  of 
Mass  Spectrometry  and  Ion  Physics,  19^,  185-200,  1976. 

7B.  H.  Mahan  and  J.  C.  Person,  "Gaseous  Ion  Recombination  Rates," 

J.  Chem.  Phys.  40,  392-401,  1964. 

8 

R.  E.  Olson,  "Absorbing  Sphere  Model  for  Calculating  Ion-Ion  Recombination 
Total  Cross  Sections,"  J,  Chem.  Phys.  56^  2979-2984,  1972. 

9 ♦ 

M.  N.  Hirsh  and  P.  N.  Eisner,  "Two-Body  Recombination  of  0,  and  02 

in  Low-Pressure  Oxygen,"  Bull.  Am.  Phys.  Soc.  17,  395,  1972. 

*®M.  N.'  Hirsh,  "Ion  Chemistry  in  Electron-Inactivated  Air:  Comparison 

of  Experiment  with  Theory,"  Bull.  Am.  Phys.  Soc.  21,  165,  1976. 
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TABLE  1.  VALUES  OF  TNO-BODY  RECOMBINATION  COEFFICIENTS  OF  ATMOSPHERIC  INTEREST 

System 

oj(300°K) 

-“-7 — 3~ 
(xlO  W/s) 

Ref 

Remarks 

1 . NO*  + N02* 

5.1  ± 1.5 

4 

mb 

1.75  ± 0.6 

5 

sa 

3.5 

A,  average  of  Ref's  445 

0.64  t 0.07 

6 

Lpfa 

2.1  ± 0.6 

7 

Sa,  value  refers  to  N0+ ( NO)  & NO^  /N02  , Ref.  3 

1.2  ± 0.3 

8 

Cal 

2.  NO+  ♦ NOj" 

8.1  ± 2.3 

4 

mb 

0.34  ± 0.12 

5 

Value  refers  to  N0+  4 N03‘/N02  (H20),  Ref.  10 

4.0 

A,  average  of  Ref.'s  445 

0.57  ± 0.6 

6 

Lpfa 

1.1  ± 0.3 

8 

Cal 

3.  NO+  ♦ O' 

4.9  ± 2.0 

4 

mb,  A 

1.9  ± 0.6 

8 

Cal 

4.  NO*  + 0{ 

5.8  l 1.0 

4 

mb,  A 

2.4  ± 0.8 

8 

Cal 

5.  02+  + N02' 

4.1  ± 1.3 

4 

i 

mb,  A 

1.2  ± 0.3 

8 

Cal 

6.  02+  ♦ NOj" 

1.3  t 0.4 

4 

mb,  A 

1.0  ± 0.2 

B 

Cal 

7,  02+  + 0- 

1.0  ± 0.4 

4 

mb,  A 

1.0  t 0.2 

8 

Cal 

8.  02+  + 02" 

4.2  ± 1.3 

4 

mb,  A 

1.0  ± 0.1 

9 

sa 

2.4  ± 0.8 

8 

Cal 

mb  = merged  beam 

sa  « stationary  afterglow 

A * values  used 

In  AIRCHEM  code 

Lpfa  * Langmuir  probe/flowing  afterglow 

Cal  * calculation 
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The  following  illustrates  the  difficulties  involved  ij  these  experi- 
ments. Recent  measurements  1 of  H30  (H20)3  + NO^  and  H^O  (H20)3  + 

N0_  (HNO,)  yield  recombination  coefficients  of  (0.55  + 0.10)  x 10  7 cm3/s 
5 * -73. 

and  (0.57  ± 0.10)  x 10  cm  /s,  respectively.  However,  the  negative 

ion  concentrations  were  mixed  with  unknown  concentrations  of: 

N02~/N02”(H20)/N03"(H20)  ions  and  NO^  (H20)/N03~ (HNOj) 2 ions,  respectively. 

The  authors  argue  that  because  the  lower  earth's  atmosphere  ought  to 
contain  mixtures  of  similar  ions,  these  values  of  a.  should  be  represen- 
tative of  the  elective  cr  in  the  earth's  stratosphere.  The  work  of 
Goldberg  et  al.  tends  to  suppgrt  £his  notion.  They  require  an  effec- 
tive a.  of  the  order  of  5 x 10  ° cni  /s  at  60  km  to  explain  ion-pair 
production  by  the  x-ray  star  Sco  X-l. 

Appeal  to  more  direct  atmospheric  measurements ^heds  little  light 
on  reducing  the  range  of  values.  Ulwick's  analysis1-5  of  the  1969  solar 
proton  event  shows  that  the  derived-effective  recombination  coefficient 
can  vary  grom^a  high  of  ~ 1.5  x 10-/  cm7s  at  45  km  to  a low  of 
~ 4 x 10  cm  /s  at  70  km.  Moreover,  he  has  estimated  that  each  of  his 
data  points  is  accurate  to  a factor  of  two,  leading  again  to  a total 
spread  of  about  a factor  of  ten. 

An  examination  into  the  sources  of  these  real  or  apparent  discrep- 
ancies is  beyond  the  scope  of  this  report.  It  is  sufficient  for  our 
purposes  to  note  that  a factor  of  10  variation  in  a nominal  value  for 
appears  to  be  reasonable. 

The  values  for  a.  used  in  the  AIRCHEM  code  are  designated  bv  an3 
"A"  in  Table  I.  All  Ather  a.  were  assigned  a value  of  2.0  x 10  cm  /s 
at  30g°K^|nd  all  the  ctj  were’* assigned  the  weak  temperature  dependence 


^D.  Smith,  N.  G.  Adams  and  M.  J.  Church,  "Mutual  Neutralization  Rates 
of  Ionospherically  Important  Ions,"  Planet.  Space  Sci.  24_,  697-703,  1976. 

12 

R.  A.  Goldberg,  W.  H.  Jones,  P.  R.  Williamson,  J.  R.  Barcus  and  L.  C. 
Hale,  "Equatorial  X-Rays  and  Their  Effect  on  the  Lower  Mesosphere," 
to  be  published  J.  Atmos.  Terr.  Phys. 

13 

J.  C.  Ulwick,  "Effective  Recombination  Coefficients  and  Lumped  Parameters 
in  the  D-Region  During  Solar  Particle  Events,"  Proceedings  of  COSPAR  Sym- 
posium on  Solar  Particle  Event  of  November  1969,  AFCRL  Special  Report 
No.  144,  571-587,  1972.  Ed.  by  J.  C.  Ulwick. 

14DNA  Reaction  Rate  Handbook,  2nd  Ed.,  March  72,  Table  241,  V. 

(Revision  No.  3,  September  1973). 
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III.  RESULTS  AND  DISCUSSION 


2 

The  AIRCHEM  code,  tailored  for  nuclear  atmospheric  disturbances, 
was  used  for  the  sensitivity  tests.  This  code  consists  of  electrons, 

15  negative  ion  species,  27  positive  ion  species  and  21  neutral  species. 

Only  the  first  hydrate  of  O2  , CO^*,  C04  , N02  and  NO^  are  considered. 

The  heaviest  and  "terminal"  negative  ion  used  in  this  code  is  N0_  (HNO,) . 
Nighttime  mass  spectrometer  flights  through  the  Dj^eg^on  have  indicated 
that  other  negative  ions  exist  in  the  atmosphere.  * Such  ions  have 
not  been  considered  in  this  model. 

A nominal  set  of  495  reactions  is  used  to  describe  the  ion  and  neutral 
chemistry.  Negative  ion  photodestruction  ha|gbeei)  included  by  considering 
the  cross  section  for  the  reaction  to  be  10  cm  for  wavelengths 
shorter  than  the  negative  ion's  electron  affinity  and  zero  elsewhere. 

Our  wavelength  region  of  interest  is  760.0  > X > 100.0  nm  and  the  solar 
flux  used  corresponded  to  midlatitude  overhead  sun  conditions.  (Recent 
computations  using  measured  cross  sections  have  not  revealed  any  sub- 
stantial changes.)  Positive  ion  photodestruction  is  not  considered  in 
this  set  of  reactions. 

The  excitation  conditions  chosen  and  fixed  for  these  sensitivity 
tests  are  thought  to  be  typical  rat^r  than  extrem^.  To  wit,  we  have 
used  for  prompt  ionization,  N = 10_^  Jjlectrons/cni  . The  delayed  ioni- 
zation is  given  by  Q(t)  = Q + t)~  ' , w^ere  t is  the  time  in  seconds 
and  we  have  selected  (^  = 10  iyg -pairs/cm  /s.  Partitioning  of  the 
charge  closely  follows  Gilmore,  except  that  N2(A3E)  is  not  considered, 
and  that  production  of  C^^A)  and  0 are  taken  as  0.25  and  1.28  per  ion- 
pair  per  second,  respectively.  Daytime  neutral  densities  are  taken  from 
CIRA  1972,  U.S.  Standard  ^mosphere  1962  and  other  literature 
sources  and  extrapolations. 

^R.  S.  Narcisi,  A.  D.  Bailey,  L.  Della  Lucca,  C.  Sherman  and  D M. 

Thomas,  "Mass  Spectrometric  Measurements  of  Negative  Ions  in  the 
D-  and  lower  E-Regions,"  J.  Atmos.  Terr.  Phys.  33,  1147-1159,  1971. 

*6F.  Arnold,  J.  Kissel,  D.  Krankowsky,  H.  Weider  and  J.  Zahringer, 

"Negative  Ions  in  the  Lower  Ionosphere:  A Mass-Spectrometric 

Measurement,"  J.  Atmos.  Terr.  Phys.  33^,  1169-1175,  1971. 

^F.  Arnold  and  D.  Krankowsky,  "Negative  Ions  in  the  Lower  Ionosphere: 

A Comparison  of  a Model  Computation  and  a Mass-Spectrometer  Measure- 
ment," J.  Atmos.  Terr.  Phys.  13,  1693-1702,  1971. 

18 

F.  Gilmore  as  quoted  by  B.  F.  Myers  and  M.  R.  Schoonover,  "Electron 
Energy  Degradation  in  the  Atmosphere:  Consequent  Species  and  Energy 

Densities,  Electron  Flux,  and  Radiation  Spectra,"  DNA  3513T,  3 Jan  75,  Table  6. 

19 

K.  S.  W.  Champion  and  R.  A.  Schweinfurth,  "A  New  Mean  Reference 
Atmosphere  for  25  to  500  km,"  AFCRL-72-0579,  2 Oct  72;  The  Mean 
COSPAR  International  Reference  Atmosphere  1972"  in  COSPAR  Inter- 
national  Reference  Atmospheres  1972,  Akademie  Verlag,  Berlin,  1972. 

20 

United  States  Committee  on  Extension  to  the  Standard  Atmosphere, 

U.S.  Standard  Atmosphere,  1962,  U.S.  Government  Printing  Office, 

Washington,  D.C.,  Dec  62. 

21 

F.  E.  Niles,  private  communication. 
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We  caution  that  the  AIRCHEM  code  has  been  validated  only  for  quiet 
or  undisturbed  daytime  conditions  and  then  only  fo^^lectron,  total 
positive  ion  and  total  negative  ion  concentrations.  This  code  has  not 
been  validated  either  at  the  level  of  excitation  or  disturbing  conditions 
or  at  the  scale  of  individual  ion  species  concentrations  used  here. 

In  the  code  each  of  15  positive  ions  is  permitted  to  recombine  with 
14  different  negative  ions,  leading  to  210  recombination  coefficients. 
Resources  did  not  permit  the  systematic  variation  of  all  possible  combi- 
nations, nor  was  this  approach  deemed  necessary.  Since  almost  all  values 
of  a.  are  identical  we  can  identify  the  dominant  contributions  to  the 
sum  defining  ou  by  noting  the  computed  dominant  ion  species  concentrations. 

For  a given  altitude  the  AIRCHEM  code  computes  ion  concentrations 
as  a function  of  tjyne  after  thg  initiation  of  the  disturbance.  At  each 
decade  of  time  (10  through  10  seconds,  inclusive)  the  dominant  positive 
ion(s)  and  dominant  negative  ions  are  identified.  Pairing  of  these 
positive  and  negative  ionic  species  then  determines  the  ion-ion  recombi- 
nation reactions  whose  coefficients  are  to  be  varied.  Table  2 shows 
these  recombination  reactions  as  a function  of  altitude.  The  corre- 
sponding rate  coefficients  were  varied  by  multiplying  a ■ by  10  and  1/10. 
(The  discyss^on  in  section  II  indicates  that  for  nominal  values  of 
a.  ~ 10  cm  /s,  decreasing  the  value  is  more  realistic  than  increasing 

tne  value.  The  increased  value  has  been  retained  for  completeness.) 

The  code  was  re-run  for  each  of  these  two  cases  and  the  ion  concentrations 
again  examined.  These  results  are  shown  in  Table  3. 

The  information  content  of  this  table  is  quite  high  and  requires  a 
bit  of  explanation.  The  altitude  is  listed  vertically  to  the  extreme 
left  and  the  decades  of  time  after  the  initiation  of  the  disturbance  are 
listed  horizontally  across  the  top  of  the  table.  The  area  within 
the  table  located  by  an  altitude  and  a time  coordinate  is  occupied  by 
two  columns  of  ionic  species,  one  positive,  the  other  negative.  To  the 
immediate  right  of  each  ion  species  is  a number  with  a superscript  and 
a subscript.  This  number  shows  the  percentage  contribution  of  that  ionic 
species  to  the  total  ion  population  for  nominal  values  of  a..  The  super- 
script (subscript)  when  added  to  the  physically  larger  number  reveals 
the  percentage  contribution  of  that  ionic  species  to  the  total  ion  popu- 
lation for  the  selected  values  of  a.  (see  Table  2)  multiplied  by 
10  (1/10).  ] 

2 

For  example,  at  an  altitude  of  60  km  and  at  10  seconds,  the  species 
NO,  constitutes  34%  of  the  total  negative  ion  population  for  nominal 
values  of  a..  For  values  of  a.  x 10  (x  1/10)  NO^  constitutes  19% 

(38%)  of  the  total  negative  ioA  population. 


22 

J.  M.  Heimerl  and  F.  E.  Niles,  "Modeling  of  Charged  Particle  Chemistry 
in  the  Stratosphere  and  Mesosphere,"  Trans.  Am.  Geophys.  Union  57,  303, 
1976. 
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TABLE  2 

ALT  (km; 
80 

70 

60 

50 

40 

30 


DOMINANT  RECOMBINATION  REACTIONS  TAKEN  FROM  BENCHMARK-76 
RESULTS;  NOMINAL  VALUES  OF  a.  HAVt  BEEN  USED. 


Dominant  Recombination  Reactions 


NO+  + 0 2~ 

NO+  + N02‘ 

N0+  + 02' 

N0+  + N03" 

N0+(H20)  + no3' 

no+(h20)2  + no3‘ 

NO+  + 02' 

NO+  + C03' 

no+(h2o)2  + no3' 
no+(h2o)3  + no3' 

N0+(H20)  + C03' 
H+(H20)3  + C03' 

H+(H2°)3  ♦ no3- 

h+(h20)4  ♦ C03~ 
H+(H20)4  + N03'(HN03) 
h+(h20)4  + N03‘(H20) 

h+ch20)4  ♦ N03"(H20) 
h+(h20)4  ♦ N03'(HN03) 
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ALT/TIME  15—i 


80  km 

°2* 

6 

72 

"14 

N0+ 

763 

3 

°2* 

9 

72 

"12 

N0+ 

0 

76 

0 

O' 

3 

16 

“3 

°2+ 

1 

20 

l 

O' 

2 

18 

“3 

°2+ 

2 

22 

2 

N02‘ 

-11 

13 

16 

N02' 

.8 

10 

12 

70  km 

V 

.1 

93 

0 

N0+ 

l 

78 

0 

°2* 

1 

89 

0 

N0+ 

0 

75 

0 

°2+ 

0 

18 

0 

°2+ 

0 

20 

0 

60  km 

V 

2 

73 

1 

N0+ 

.3 

71 

0 

V 

1 

56 

1 

N0+ 

_4 

56 

0 

co3- 

.3 

15 

1 

°2+ 

1 

11 

0 

co3- 

.2 

30 

0 

no+(h2o) 

_1 

14 

0 

°2+ 

1 

9 

0 

50  km 

co3- 

.2 

42 

1 

no+(h2o) 

_7 

29 

0 

C03- 

2 

55 

1 

no+(h2o) 

.9 

36 

0 

o2- 

1 

19 

0 

N0+ 

3 

27 

0 

no3- 

4 

15“ 

2 

h+(h20)3 

.1 

18 

.1 

no3-* 

2 

12 

0 

h+(h2o)2 

2 

11 

0 

no*(h2o)2 

.3 

14 

1 

co4' 

0 

10 

0 

h+(h20)2 

5 

11 

"1 

NO* 

4 

9 

•1 

40  km 

co3- 

.3 

47 

“3 

h+(h2o)4 

.28 

58 

7 

C03'(H20) 

-15 

31 

9 

h+(h20)4 

.14 

67 

4 

C03'(H20) 

-15 

31 

8 

h+(h2o)3 

.6 

18 

1 

co3- 

4 

31“ 

0 

h+(h20)3 

-3 

19 

0 

no3'(h2o) 

9 

11 

•4 

N0+(H20)2 

19 

14 

•7 

no3'(h2o) 

18 

27 

"6 

no+(h2o)2 

9 

8 

*3 

30  km 

C03'(H20) 

-31 

69 

11 

h+(h20)4 

2 

71 

4 

co3'(h2o) 

-19 

58 

5 

h+(h20)4 

1 

75 

2 

N03'(H20) 

28 

19 

•7 

h+(h20)5 

0 

23 

0 

no3'(h2o) 

19 

29 

•7 

h+(h20)5 

0 

23 

0 

no3‘(hno3) 

_5 

10 

2 

TABLE  3.  PERCENTAGE  CONTRIBUTION  Of  LISTEO  ION  TO  THE  TOTAL  ION  POPULATION  AS  A FUNCTION 
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’02s 

103s 

1G4  s 

V 

15 

52 

”19 

N0+ 

0 

91 

0 

no2" 

-29 

41 

13 

N0+ 

0 

95 

0 

no2' 

-26 

44 

6 

N0+ 

0 

88 

0 

no2' 

-20 

24 

27 

V 

25 

32 

"li 

V 

27 

30 

"7 

o' 

6 

22 

“0 

0* 

12 

15 

"5 

N03' 

18" 10 

3 

N03- 

.7 

10 

4 

o' 

7 

8 

"2 

V 

3 

68 

"9 

N0+ 

0 

78 

0 

V 

16 

42 

"17 

N0+ 

0 

59 

0 

no3‘ 

-34 

57 

10 

N0+ 

0 

31 

0 

• 

on 

O 

u 

4 

18 

•7 

no+(h2o) 

0 

10 

0 

N0j* 

-20 

34 

25 

no+(h2o) 

0 

25 

0 

°2* 

14 

19 

"5 

N0+(H20)2 

0 

25 

1 

N03' 

1 

5 

17 

co3' 

5 

12 

"5 

no2' 

11 

15 

"4 

no+(h2o) 

0 

23 

0 

co3- 

8 

6 

0 

c03' 

10 

41 

0 

N0+ 

32* 1 

0 

N03* 

-26 

63 

9 

no+(h2o)2 

.4 

38 

l 

N03" 

-24 

71 

6 

h+(h20)3 

- 3 

38 

2 

no3' 

-15 

34 

4 

N0+(H20) 

.2 

23 

0 

co3- 

19 

16 

"4 

no+(h2o) 

0 

18 

0 

no2" 

10 

9 

"3 

h+(h2o)4 

.1 

25 

1 

V 

5 

15 

0 

mo+(h20)2 

18" 1 

0 

h+(h2o)3 

-1 

14 

0 

C03- 

7 

6 

1 

no+(h2o)2 

0 

24 

-2 

H+(H20)2 

-l 

14 

0 

h+(h20)2 

3 

9 

0 

no+(h2o) 

1 

9 

"1 

co3- 

4 

41 

“4 

h+(h20)3 

-13 

43 

8 

N03- 

.6 

47 

3 

h+(h20)3 

- -22 
53 

9 

no3- 

_ 4 

49 

1 

h+(h2o)3 

-25 

62 

6 

no3- 

-10 

37 

5 

N0+(H20) 

2 

32 

“8 

N03'(H20) 

2 

32 

"l 

N0+(H20) 

26 1 3 
"7 

no3'(h2o) 

0 

43 

"2 

no+(h2o) 

17 

20 

"6 

no3- 

(H90)  171 
c “l 

no+(h2o)2 

2 

13 

“3 

C03- 

7 

15 

"2 

no+(h2o)2 

6 

11 

"3 

no+(h2o)2 

8 

8 

"1 

N03‘(H20) 

h+(h20)4 

.5 

77 

”4 

N03'(H20)  7l" 

M^(Ho0).  76" 

^ 4 0 

C03'(H20)  12j6 

h+(h20)3 

.1 

21 

”1 

N03'(HN03)  14’ 7 

H^(H,0),  21° 

3 0 

N°3" 

12’1 

"l 

N0.,"(HN0,)  38’ 16 
J 3 6 

h+(h20)4 

-3 

7*.3 

C03'(H20) 

2 

34 

"5 

h+(h2o)5 

0 

24 

0 

no3'(h2o) 

12 

30 

"5 

OF  ALTITUDE , TIME  AND  CHANGES  IN  SELECTED  VALUES  OF  . (See  text  for  discussion.) 


IS 


The  values  in  Table  3 are  given  in  percent  of  total  ion  population 
to  account  for  the  fact  that  changing  values  of  a.  can  alter  the  total 
ion  concentrations.  Thus  changes  in  the  species  Attributions  are  more 
easily  read  from  this  table.  Code  predicted  absolute  values  can  be 
recovered  for  each  entry  in  Tabl^  3 by  multiplying  by  the  appropriate 
total  ion  concentrations  (in  cm  listed  in  Table  4.  The  key  "MAX," 
"NOM,"  and  "MIN"  correspond  to  the  cases  a.  x 10,  4.  and  a.  x 1/10, 
respectively.  ^ ^ 


This  version  of  the  AIRCHEM  code  computes  a.'  = a.  + a,_  x [M] , 
where  a.  is  the  two  body  ion-ion  recombination  cAefficlent  (of  interest 
in  this J study),  a,„  is  the  three  body  ion-ion  recombination  coefficient 
and  [M]  is  the  total  neutral  concentration  at  a given  altitude.  Since 
only  a.  has  been  varie^j  the  effect  of  a_g  x [M]  must  be  accounted  for. 
From  Cole  and  Pierce's  "selected  value"  for  theirfiiog-ion  recombina- 
tion coefficient  at  ground  level  the  value  8 x 10''cm  /s  can  be  derived 
for  a,_.  This  value  is  such  that  a.'  » a.  except  at  30  km  where 
a./10  * a,_  x [M] . This  case  was  re-compAted  with  a,B  set  equal  to  zero. 
Trivial  changes  (0-3  percentage  points)  were  found  inBthe  computed  ion 
distributions  of  Table  3 and  small  changes  (12-15%)  were  found  in  the 
total  ion  concentrations  of  Table  4.  Thus,  no  corrections  for  the 
inclusion  of  the  three  body  ion-ion  recombination  have  been  made  to  the 
values  listed  in  Tables  3 and  4. 


In  Table  4 at  70  and  80  km  altitudes,  large  differences  exist 
between  the  sum  of  the  positive  and  sum  of  the  negative  ions.  Since 
charge  is  conserved  in  the  AIRCHEM  code,  this  difference  defines  the 
electron  density,  and  shows  that  electron-ion  recombination  does  dominate 
at  the  higher  altitudes.  As  we  progress  lower  in  altitude  (and/or  later 
in  time)  the  sum  of  the  negative  ions  approaches,  then  approximately 
equals  the  sum  of  the  positive  ions.  Thus  for  our  choice  of  excitation 
conditions,  ion-ion  recombination  becomes  important  at  60  km  and  dominates 
below  50  km. 


Referring  to  Table  3 we  find  no  case  which  shows  more  than  a 
factor  of  about  two  change  in  a dominant  ion  species.  Specifically  see 
H (l^O)^  at  40  kip  and  CO^  (^0)  at  30  km  both  at  10  seconds.  These 
changes  are  for  the  less  likely  increase  in  a.;  the  more  likely  decrease 
indicates  rather  small  changes  in  the  ionic  distributions  for  altitudes 
60  km  and  below. 


The  AIRCHEM  code  has  been  used  as  a tool  to  study  the  sensitivity 
of  distributions  to  changes  in  selected  values  of  individual  two  body  ion- 
ion  recombination  coefficients.  We  caution  that  the  code  has  not  been 
validated  for  the  conditions  or  in  the  detail  used  herein,  and  that  the 
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R.  K.  Cole,  Jr.  and  E.  T.  Pierce,  "Electrification  in  the  Earth's 
Atmosphere  cor  Altitudes  Between  0 and  100  Kilometers,"  J.  Geophys. 
Res.  70,  2735-2749,  1965. 
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Total  Total  Total  Total  Total  Total  Total  Total  Total  Total 


TABLE  4.  COMPUTED  TOTAL  ION  POPULATION  CONCENTRATIONS  (in  cm"J)  AS  A FUNCTION  OF  ALTITUDE,  TIME  AND  CHANGES 
IN  SELECTED  VALUES  OF  . (See  text  for  discussion) 


ion  species  used  are  a sub-set  of  those  detected  in  the  atmosphere. 

With  these  points  in  mind  we  find  the  results  of  this  limited  study 
indicate  that  for  altitudes  lower  than  60  km,  one  extreme  value  of  the 
ion-ion  recombination  coefficient,  o.  x 10,  can  change  ionic  populations 
by  as  much  as  a factor  of  two,  in  isolated  cases.  The  other,  more 
probable  extreme,  a.  x 1/10,  causes  much  smaller  changes  in  the  computed 
ion  populations.  ^ 
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